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ABSTRACT 

We calculate simultaneously the radio and optical luminosity evolutions of quasars, and the distri- 
bution in radio loudness R defined as the ratio of radio and optical luminosities, using a fiux limited 
data set containing 636 quasars with radio and optical fiuxes from White et al. We first note that 
when dealing with multivariate data it is imperative to first determine the true correlations among 
the variables, not those introduced by the observational selection effects, before obtaining the indi- 
vidual distributions of the variables. We use the methods developed by Efron and Petrosian which 
are designed to obtain unbiased correlations, distributions, and evolution with redshift from a data 
set truncated due to observational biases. It is found that the population of quasars exhibits strong 
positive correlation between the radio and optical luminosities. With this correlation, whether in- 
trinsic or observationally induced accounted for, we find that there is a strong luminosity evolution 
with redshift in both wavebands, with significantly higher radio than optical evolution. We conclude 
that the luminosity evolution obtained by arbitrarily separating the sources into radio loud {R > 10) 
and radio quiet {R < 10) populations introduces significant biases that skew the result considerably. 
We also construct the local radio and optical luminosity functions and the density evolution. Finally, 
we consider the distribution of the radio loudness parameter R obtained from careful treatment of 
the selection effects and luminosity evolutions with that obtained from the raw data without such 
considerations. We find a significant difference between the two distributions and no clear sign of 
bi-modality in the true distribution for the range of R values considered. Our results indicate there- 
fore, somewhat surprisingly, that there is no critical switch in the efficiency of the production of disk 
outflows/jets between very radio quiet and very radio loud quasars, but rather a smooth transition. 
Also, this efficiency seems higher for the high-redshift and more luminous sources in the considered 
sample. 

Subject headings: quasars: general - methods: data analysis - methods: statistical - galaxies: active - 
galaxies: jets 



1. INTRODUCTION 

The optical emission of quasars, or active galactic nu- 
clei (AGN) is dominated by the radiation of the plasma 
accreting onto super massive black holes, while the radio 
emission is dominated by the plasma outffowing from 
the black hole/ accretion disk systems. Hence different 
but complementary information can be gathered in both 
photon energy ranges regarding the cosmological evolu- 
tion of AGNs and its relation to structure formation in 
the Universe. It is therefore important to analyze in de- 
tail redshift distributions of quasars in both frequency 
regimes, investigating carefully any possible differences 
between these two. 

The rapid evolution of active galaxies identified in ra- 
dio catalogs as 'quasi-stellar radio sources' (or QSRs) 
in the redshift range z < 2 was established soon af- 
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ter their discovery (e.g. iSchmidt 1 119671 ). Subsequent 
optical discoveries of similar sources, most of which 
had no detectable radio emission, lead to the emer- 
gence of the class of 'radio quiet quasars ' (or 'quasi- 
stellar objects'; QSOs for short e.g. lOsterbrock I 
(1989)). These optically-selected sources also showed 
similar strong evolutionary trends, similar to the radio- 
selected ones. These evolutions are modeled as den- 
sity evolution, luminosity evolution or a combination 
of the two in numer ous wor ks (e.g. iSchmidtl 119671: 
Petrosian 1973; Marsh all et al.ll l983: Du nlop fc Peacocp 
1990; Malo nev fc PetrosianI 11999 : Willott et al.l I2QQ1D 
and can be designated as the evolution of the luminosity 
function (LF, for short). 

By now the evolution of the LF has been described 
not only for optical and radio luminosities but also 
for X-ray, i nfrare d, and b olomet r ic luminosities (e.g. 
Ueda et all I2QQ3I: Richard s et al.l 120061: iMatute et al.l 
20061 : iHopkins et al . 2007; Croom et al.ll2009f ). Most of 



these studies have treated the evolution with a bi-variate 
function ^i(L^, z), where Li is the luminosity (or, in this 
case, luminosity spectral density) in some photon energy 
range, e.g. Li = Lopt or Lrad- The shape of the LF 
and its evolution are usually obtained from a flux lim- 
ited sample fi > fm,i with Li = Att d\{z) {1/ Ki{z)) fi, 
where is the luminosity distance and Ki{z) stands for 
the K-correction. For a power law emission spectrum of 
index £i defined as /i oc one has i^i(z) = {l-\-zy~^\ 
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However, because no matter how a quasar is discov- 
ered, optical observations are required for determination 
of redshift, then the flux limit of optical observations 
{fm,opt) and the optical luminosity enter the picture, so 
that one now must consider the joint LF and its evolu- 
tion ^(Lopt, i^i, a tri-variate function with Li = L^ad 
or Lx, for example. 

In general, the first step required for investigation of a 
multivariate distribution is the determination of whether 
the variables of the distributions are correlated or are 
statistically independent. For example, in the case of 
a single LF the correlation between L and z is what we 
call luminosity evolution, and independence of these vari- 
ables would imply absence of such evolution. Mathemat- 
ically, independence means that the function is separable 
= X p{z)^ in which case one is left with 

the determination of a single variable LF ijjiiLj) and the 
density evolution p{z). As shown by iPetrosian I (11992^ 
the most exact nonparametric method for this task from 
a fl ux limited (or a mo re generally truncated) sample is 
the iLvnden-Belll (|1971[ ) method. However, this simple 
and elegant method cannot be used for cases when vari- 
ables are correlated ( e.g. when there is luminosity evo- 
lution). lEfron fc PetrosianI (|1992L 11999 ) (EP for short) 
developed new methods for determination of the exis- 
tence of correlation or independence of the variables from 
a flux limited and more generally truncated data set, 
and prescribed how to remove the correlation by defin- 
ing new and independent variables [say L[ = Li/gi{z) 
and z, where the function gi{z) describes the luminosity 
evolution] and then how to determine the mono-variate 
functions il^i{L'-) and p{z). Thus, one can determine both 
the luminosity and density evolutions gi{z) and p{z)^ as 
well as the LF at any redshift 

In the case of quasars with the optical and some other 
band luminosity, we have at least a tri-variate function. 
In this case one must determine not only the correla- 
tions between the redshift and individual luminosities 
(i.e. the two luminosity evolutions) but also the possible 
correlation between the two luminosities, before individ- 
ual distributions can be determined. Knowledge of these 
correlations and distributions are essential for not only 
constraining robustly the cosmological evolution of ac- 
tive galaxies, but also for interpretation of related obser- 
vatio ns, such as the ex trags:alactic background radiation 
(e.g. Singal e t al.]l2QlQl : IHopkins et al.ll2QlQl) . 

Another related aspect of this subject, which has at- 
tracted considerable attention over the years, is the dis- 
tribution of the'radio-loudness parameter' for the quasar 
population, and the distinction between so-called 'radio 
loud' (RL for short ) and 'radio quiet' (RQ for short) 
quasars. The question of whether there are two dis- 
tinct populations was addressed soon after the discov- 
ery of quasars using small samples with radio flux lim- 
its greater than one Jy. Initially it was found that the 
distribution of the ratio of radio to optical luminosity 
R = ^rad/^opt, the SO Called "radio loudness parame- 
ter", was a fairly broad power law with index j3R ~ —2.3 

^ It should be noted that here we assume that the shape of 
the LF is constant; e.g. power law indices describing the LF are 
independent oi z. In general, shape variations can affect the test 
of independence. For a sufficiently large sample the importance of 
these effects can be determined and accounted for. This is beyond 
the scope of this paper. 



in th e range 2.8 < \ogR < 5.2 (iSchmidt II1972I : IPetrosianI 
Il973f ). At that time this ratio was defined for the ra- 
dio luminosity at z^rad = 0-5 GHz and optical luminos- 
ity at 2500 A (or the frequency Uopt = 1.2 x 10^^ Hz). 
Nowadays it is defined with radio luminosity at 5 GHz 
so that for a mean radio spectral index ^rad ^ 0.6 the 
old data would be in the range 2.2 to 4.6 of the mod- 
ern definition of log i^H Later, however, the survey lim- 
its were extended to much lower fluxes (specially in the 
radio domain), and this has resulted in a much wider 
range of the ratio that extends to values well below one, 
namely — 3 < \ogR < 5. Within this bro ader rang:e, weak 
hints of the bi-modality described by (jKellerman et al.l 
Il989l ) suggested that logR = 1 could be chosen as the 
radio loud/quiet demarcation value. Using this value 
for the division between RL and RQ quasars, the dif- 
ferences between the two classes have been investigated, 
including the possibility of distinct cosmolo gical evolu- 
tion of the RL and RL populations (e g Mille r et al.lll990l : 
iGoldschmidt et al.l 119991 : iJiang et al. 2007). Stih, the 
more recent analyses of different samples of objects re- 
ported in the literature so far gave rather inconclusive 
results on whether any bi-modiality in the distribution 
of the radio loudness parameter for quasars is inherent 
in the population (see llvezic et al.|[2002l : iCirasuolo et al.l 
[20031 : llvezic d~ZI[200l 

There have been many papers dealing with this ratio 
and RL vs RQ issue, as well as luminosity ratios at other 
wavelengths, e.g. IR/radio, Optical/X-ray etc. However, 
none of these works have dealt with the intrinsic distri- 
bution (and/or evolution) of the ratio, which is related 
to the tri-variate LF ^(Lopt, ^rad, ^) b}|3 

Gr{R, z) = ^(i^opt, R ^opt, z) Lopt dLopt 
Jo 

These works did not take the observational selection ef- 
fects properly into consideration, nor did they address 
the correlations between the radio and optical luminosi- 
ties. Neglecting these effects when attempting to deter- 
mine the distribution of radio loudness is usually given 
the justification that the ratio is essentially independent 
of cosmological model and redshift (as long as the K- 
corrections are the same). The broad distribution of ob- 
served ratios obtained in this way (see Figure [14] below) 
deviates from a simple power law and may even have a 
hint of bi-modality, seemingly justifying at face value the 
choice of log ~ 1 as the separation point between RL 
and RQ sources. However, as shown in Appendix A even 
in the simplest cases the observed distribution (and its 
moments) could be very different from the intrinsic ones. 

^ The fiducial cosmological model used at that time, namely 
the Einstein-De Sitter model, was also different than the currently 
accelerating models. However, this will affect the values of the 
individual luminosities but not the ratio R. We also note that 
some other authors use B-band optical fluxes in defining the radio 
loudness parameter, but the difference is not large, resulting in a 
change in the ratio by a factor of 1.33, for the assumed optical 
spectral index Sopt = —0.5. 

^ Equation [T] arises because by definition f GR(R,z)dR = 
f f ^(Lopt, L Y z) dL opt dL^Q^fi, and from the definition of i?, 
c^^rad = ^opt dR and dLopt = -(Lrad/^^) dR- 
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Thus, for determination of the true distribution of R the 
data truncations must be determined and the correla- 
tions between all variables must be properly evaluated. 

Our aim in this paper is to take all these effects into 
account in determination of the evolution of optical and 
radio luminosities and their ratio and to find their dis- 
tributions. In ^ we describe the data we use. In ^ 
we provide an overview of the procedure used. In ^ 
we present our results on the correlations and evolu- 
tions of the LFs. In ^ we describe the density evolu- 
tion and the luminosity density evolutions, while in ^ 
we calculate the LF corrected for luminosity evolution, 
which we call the "local" LF. Finally, in ^ we evalu- 
ate the distribution of radio loudness, R. This work as- 
sumes the standard ACDM cosmology throughout, with 
Ho 71kms-^Mpc-\ 0.7 and Qrn = 0.3. 

2. DATASET 

In order to evaluate the luminosity evolution in both 
radio and optical, and to separate and compare these 
effects, we require a data set that has both radio and 
optical fluxes to reasonable limits and across a range 
of redshifts, that contains a significant number of both 
RL and RQ objects. The overlap of the FIRST bright 
quasar radio survey with the Automatic Plate Measuring 
Facility catalog of the Pal omar Observatory S ky Survey 
(POSS-I), as presented bv lWhite etHI ([2000'). is such a 
data set. It contains 636 objects with optical R band 
optical magnitudes, 1.4 GHz total and peak pixel fiuxes, 
and spectroscopic redshifts. The survey has a limiting 
R band magnitude of 17.8 or fm,R = 0.22 mJy, a limit- 
ing peak pixel 1.4 GHz fiux of 1 mJy, and redshifts that 
range from 0.02 to 3.425. Figures [1] and [2] show the radio 
and optical luminosities versus redshifts of the quasars 
in the survey, assuming the standard K-corrections for 
power laws with optical and radio indices Sopt =0.5 and 
and Srad = 0.6 respectively (where fi ex Figure [3] 

shows the radio loudness parameter R versus redshift for 
the dataset. 

This sample spans a very wide range of luminosities 
(5 dex in optical and 7 dex in radio) with a significant 
number of sources in the range 0.1 < R < 10^ (with 
336 RL and 300 RQ). Therefore, it is well suited for our 
analysis here. We have examined some other combined 
radio and optical survey data sets and found them to 
be not as well suited for this analysis. For example, the 
combined FIRST r adio survey with the 2d F optical sur- 
vey as reported in iCirasuolo et al.l (|2003[ ) features only 
12 RQ objects (of 113 total), and the combined FIRST 
wi th the Larg e Bright Quasar Survey (LBQS) as reported 
in iHewett et al.i (f2001i) has only 77 objects and differ- 
ent optical fiux limits for the various fields, making the 
method employed here cumbersome. Of course a much 
larger sample could be achieved combining optical data 
from the Sloan Di gital Sky Su r vey ( SPSS) with FIRST 
radio data, as in Uiang et al.l (120071). Forming such a 
dataset, however, would necessitate care in associating 
optical sources with radio ones, given the very large num- 
ber of optical sources. Issues such as how large of a radius 
to allow for an association and what to do about mul- 
tiple matches are important considerations without im- 
mediate answers. It is our intention to demonstrate the 
techniques employed here with a well established smaller 
dataset before moving on to a more comprehensive one 
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Fig. 1. — The 2 500 A rest frame a bsolute luminosity density for 
the quasars in the White et al.' (2000^) dataset used in this analysis. 
To obtain the 2500 A luminosity density we convert from observed 
i?-band magnitude to flux at the integrated center band frequency, 
and assume an optical spectral index of 0.5 and the luminosity 
distance obtained from the redshift with the standard cosmology 
and the standard K-correction. The crosses are the RL objects 
while the diamonds are the RQ. 
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Fig. 2. — The 1.4 GHz rest frame absolute luminosity density for 
the quasars in the White et al. ( 2000) dataset used in this analysis. 
To obtain the 1.4 GHz luminosity density we use the luminosity 
distance obtained from the redshift and the standard cosmology 
and the standard K-correction. We assume an radio spectral index 
of 0.6. The crosses are the RL objects while the diamonds are the 
RQ. 

involving SDSS optical data. 

3. GENERAL REMARKS ON CORRELATIONS IN 
LUMINOSITY FUNCTIONS 

The LF gives the number of objects per unit comov- 
ing volume V per unit source luminosity, so that the 
number density is dN/dV = / dLi'^i{Li^ z). To examine 
luminosity evolution, without loss of generality, we can 
write a LF in some waveband i as 

^i{Li,z) = f{z)tljlLi/glz),r]i)/gi{z), (2) 

where gi{z) and f{z) describe the luminosity evolution 
and comoving density evolution with redshift respec- 
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Fig. 3. — The redshift distribution of the ratio R of rest frame 
abs olute lumin osities at 5 GHz and 2500 A for the quasars in 
the lWhite et al. (2000) dataset used in this analysis. The 5 GHz 
luminosity is obtained from the 1.4 GHz luminosity assuming a 
radio spectral index of 0.6. The crosses are the RL objects while 
the diamonds are the RQ. 

tively and r]j stands for parameters that describe the 
shape (e.g. power law indices and break values) of the i 
band LF (we use the normalization tlJi{Li)dLi = 1)0 
In what follows we assume a non-evolving shape for the 
LF (i.e. = const, independent of L and z), which 
is a good approximation for determining the global evo- 
lutions. Once these are determined this hypothesis can 
be tested and results amended. However, for more com- 
plicated functional forms with variable r]j^ e.g. for lu- 
minosity dependent density evolution, the determination 
of the variations will require a large sample with signif- 
icant numbers of objects in reasonably narrow redshift 
and luminosity bins. 

Given this assumption then once the luminosity evolu- 
tion gi{z) is calculated, the density evolution f(z) and lo- 
cal LF i^i{L^) = tlJi{Li/gi{z))/gi{z) can be determined^ 
We consider this form of the LF for luminosities in dif- 
ferent bands, allowing for separate (optical and radio) 
luminosity evolution. 

1. As is often done, one might naively assume that 
the joint LF ^(Lopt, i^rad, ^) is separable into two forms 
like equation [2] with a common density evolution. How- 
ever, as discussed in ^ because the optical and ra- 
dio luminosities of the quasars are, in general, highly 
correlated, the simultaneous determination of the LFs 
of both requires care. The first step in this procedure 
should be to determine the degree and form of the cor- 
relation between the optical and radio luminosities. As 
described below, the EP method allows us to determine 
whether any pair of variables are independent or corre- 
lated. Once it is determined that they are correlated 

^ There are in principle other possible parameters, e.g. the spec- 
tral indices. We can ignore them for the purposes of the analysis 
here on the assumption that they either do not evolve strongly with 
redshift or are not strongly correlated with any of the luminosities 
in question. 

The method developed by EP that we shall use below actually 
gives the cumulative functions cr(< z) = p{z')[dV{z') / dz']dz' 
and (/)(> L') = i^(L'') dL" . The differential functions p and 
are obtained by differentiation. 



one should seek a coordinate transformation to define a 
new pair of variables which are independent. This re- 
quires a parametric form for the transformation. One 
can define a new luminosity which is a combination of the 
two; we can define a "correlation reduced radio luminos- 
ity" Lcrr = Lrg,d/F{Lopt/Lfid), whcrc the function F de- 
scribes the correlation between Lrad and Lopt and L fid is 
a fiducial luminosity taken here to be 10^^ erg sec~^ Hz~^ . 
This is a convenient choice for L fid as it is lower than the 
lowest 2500 A luminosity considered in our sample, but 
results do not depend on the particular choice of numer- 
ical value. For the correlation function we will assume a 
simple power law 



(3) 



where a is a bulk power law correlation index to be de- 
termined by a fit to the data. This is essentially a coordi- 
nate rotation in the log- log luminosity space. As shown 
in ^ below, EP also prescribe a method to determine a 
best fit value for the index a which orthogonalizes the 
new luminosities. Given the correlation function we can 
then transform the data (and its truncation) into the new 
independent pair of luminosities (I/opt and Lcrr), whose 
distribution can be represented as 



^(L 



opt 1 L/QYYl ^) 



p{z) X V^opt(^opt/^opt,^opt)/^opt 
X V^crr(^crr/^crr,^crr)/^ err- 



(4) 



2. The next step is determination of the two indepen- 
dent luminosity-redshift correlation functions ^opt and 
^crr which describe the luminosity evolutions. The pro- 
cedure for determination of these functions is similar to 
the ones for removing the correlations between the lumi- 
nosities except now we make coordinate transformations 
in the Lopt — z and Lcrr — z spaces. We assume simple 



forms 



gi{z) = il + z)' 



(5) 



SO that L[ = Li/gi{z) refer to the local {z = 0) luminosi- 
tiesQ The full procedure is detailed in §31 

3. The density evolution function p{z) is determined 
by the method shown in EP (see ^ below). Once all 
correlations are removed we end up with a local separable 
LF as in equation [H 

4. The local LFs of uncorrelated luminosities L^p^ and 
L^^^ can then be used to recover the local radio LF by a 
straight forward integration over L'^^^ and the true local 
optical LF as 



Jo 



^rad(i^rad) 



V^opt(^opt) 



LI 



rad 



dLl 



opt 



As stated above this procedure can be used for the de- 
termination of the radio LF at any redshift, from which 

This is an arbitrary choice. One can chose any other fiducial 
redshift by defining gi{z) = [(1 + z)/{l + Zfid)]^^. 
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one can deduce that the radio luminosities also undergo 
luminosity evolution with 

^rad(^) = Qcriz) X [^optl^)]"" (7) 

(cf equation [3|) 

5. Similarly we can determine the local distribution of 
the radio to optical luminosity ratio, R' = L'^^^/ L'^^^ = 

^crr ^ ^opt" ^ ^ ^fid "5 ^S 



(8) 

and its evolution 

giiz) = gUz) X [go^zT-' = ^— (9) 

^opt 

4. CORRELATION FUNCTIONS 

We now describe results obtained from the use of the 
procedures described in ^on the data described in ^ 
Here we first give a brief summary of the algebra involved 
in the EP method. We follow c losely the steps described 
in iMalo ney fc PetrosianI () 19991 ). This method uses the 
Spearman rank test to determine the best-fit values of 
parameters describing the correlation functions using the 
test statistic 

(10) 

to test the independence of two variables in a data set, 
say {xj^Hj) for j = 1, . . . ,n. Here Rj is the y rank of 
the data point j in a set associated with it. For a un- 
truncated data (i.e. data truncated parallel to the axes) 
the associated set of point j includes all of the Xk < Xj . 
If the data is truncated one must form the associated 
set consisting only of those points of lower x value that 
would have been observed if they were at the x value of 
point j given the truncation. As an example, if we have 
one sided truncations as in Figures [1] and [21 then the as- 
sociated set Aj = {k : yk > yj, y'^ < yj}^ where y^ 
is the limiting y value of data point i (see EP for a full 
discussion of this method). 

If {xj , yj ) were independent then the rank IZj should 
be distributed uniformly between and 1 with the ex- 
pectation value and variance £j = (1/2) (j + 1) and 
Vj = (1/12) (j^ + 1), respectively. Independence is re- 
jected at the n a level if \r\ > n. To find the best 
fit correlation the y data are then adjusted by defining 
yj ~ yj/-^i^j) rank test is repeated, with differ- 

ent values of parameters of the function F. 

4.1. Radio- Optical Luminosity Correlation 

The radio and optical luminosities are obtained from 
radio and optical fluxes from a two flux limited sample 
so that the data points in the two dimensional flux space 
are truncated parallel to the axes which we consider to 
be untruncated. Since the two luminosities have essen- 
tially the same relationship with their respective fluxes, 
except for a minor difference in the K-correction terms. 
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Fig. 4. — The value of the r statistic as given by equation llOlas a 
function of a for the relation Lrad oc (Lopt)"", where Lopt and Lrad 
are the optical and radio luminosities, respectively, for the quasars 
in the dataset. The 1 a range for the best fit value of a is where 
I r I < 1. It is seen that the observed optical and radio luminosities 
are strongly positively correlated, with a linear or slightly higher 
power law relation. 

we can consider the luminosity data to also be untrun- 
cated. In that case as mentioned above the determina- 
tion of the associated set is trivial and one is dealing with 
the standard Spearman rank test. Assuming the corre- 
lation function between the luminosities F{x) = we 
calculate the test statistic r as a function of a. Figured] 
shows the absolute value of the r vs a, from which we get 
the best flt value of a = 1.3 with one a range ±0.2. As 
expected a is near unity, and the value a=l is not ruled 
out with a high significance. As discussed in Appendix 
B, this correlation may be inherent in the population or 
may be an artifact of the fiux limits and wide range of 
redshifts, although this particular point is not important 
for the analysis going forward, as the rotation to Lcrr is 
a technique to achieve independent variables (Lopt and 
Lcrr) in the context of the data present to recover the 
inherent redshift evolutions. 

4.2. Luminosity- Redshift Correlations 

We now describe our results on determination of the 
luminosity evolution, i.e. the luminosity-redshift corre- 
lation functions gi{z), which according to equation [5] re- 
duces to determination of the values of the indices ki. 
The basic method for determining the best fit ki is the 
same as above but in this case the procedure is more 
complicated for several reasons. First, as evident from 
Figures [T] and [2] the Li — z data are heavily truncated 
due to the flux limits. Second, we now are dealing with 
a three dimensional distribution (Lcrr, ^opt, ^) and two 
correlation functions [^crr(^) and ^opt(^)]- 

Specifically, since we have two criteria for truncation, 
the associated set for each object includes only those ob- 
jects that are sufficiently luminous in both bands to ex- 
ceed both ffux minima for inclusion in the survey if they 
were located at the redshift of the object in question. 
Consequently, we have a two dimensional minimization 
problem, because both the optical and correlation re- 
duced radio evolution factors, ^opt(^) = (1 + z)^°p^ and 
^crr(^) = (1 + zY""" ^ come into play, as the luminosity 
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cutoff limits for a given redshift are adjusted by powers 
of /Copt and /Crad too. 

We form a test statistic Tc, 



opt 



r^^^ where 



Topt and Terr are those evaluated considering the objects' 
optical and correlation reduced radio luminosities, re- 
spectively. The favored values of /Copt and /Ccrr are those 
that simultaneously give the lowest Tcomh and, again, we 
take the la limits as those in which Tcomh < 1- For 
visualization. Figure [5] shows a surface plot of Tcomb 

We have verified this method with a simulated Monte 
Carlo data set in which objects are distributed in red- 
shift and given randomized luminosities in accordance 
with set optical and radio evolutions. The algorithm can 
recover the evolutions correctly provided that they aren't 
wildly different, i.e. one very positive and the other very 
negative. 

Figure [6] shows the best fit values of kopt and kcrr and 
taking thel and 2 a contours. Results are shown for the 
entire dataset taken as a whole, and also with the data 
split into the RL and RQ subsets. The radio luminosity 
evolution itself can be recovered by equation [71 

Given the tight constraints achieved when the dataset 
is considered as a whole, and the sharp bifurcation when 
the set is split into the RL and RQ populations, it is 
evident that splitting the population before determina- 
tion of the luminosity evolutions introduces a bias into 
the determinations. This is expected because differing 
evolutions will have a strong effect on the likelihood that 
objects at a given redshift will fall on the RL or RQ side 
of the division according to the standard definition used 
here, and the data for each set will be artificially trun- 
cated along = 10 as a function of the evolutions. 

We see that positive evolution in both radio and opti- 
cal wavebands is favored. The minimum value of Tcomb 
favors an optical evolution of /Copt = 3.0 and a radio evo- 
lution of /crad = 5.4, but uncertainty at the la level allows 
the range of /Copt from 2.5 to 3.25, and /crad from 5.3 to 
5.75. Therefore, we conclude that quasars have under- 
gone a significantly greater radio evolution relative to 
optical evolution with redshift. In the above analysis we 
have assumed sharp truncation boundaries and that the 
data is complete above the boundaries. As discussed in 
^this may not be the case for the FIRST radio data. If 
an estimate of the uncertainty from the consideration of 
possible radio incompleteness at faint fluxes is included, 
the favored optical range enlarges to kopt from 1.25 to 
3.75, with a slightly lower best fit value of kopt = 2.0. 
Due to combined effects on /copt and /Ccrr, the value of k^ad 
is not much affected by allowing for possible radio incom- 
pleteness at faint fiuxes, perhaps counter intuitively. 

5. DENSITY EVOLUTION 

Next we determine the density evolution f(z). One can 
define the cumulative density function 



aiz)= f dz) 
Jo 



dz 



(11) 



which , fohowing iPetrosianl ()1992h based on iLynden-Belll 

can be calculated by 




Fig. 5. — Surface plot of the value of Tcomb for the dataset as 
a whole showing the location of the minimum region where the 
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Fig. 6. — The la and 2a contours for the simultaneous values of 
fcopt and fccrr where the optical and correlation reduced radio lumi- 
nosity evolutions are gopt(z) = (1 + z)^°p^ and gcrr(z) = (1 + z)^"^" . 
The best fit radio luminosity evolution can be reconstructed from 
9rad = 9cTT X Rcsults are shown for the data set evaluated 

as a whole (solid contours), and for the RL (dash-dot contours) 
and RQ (dashed contours) populations evaluated separately. It is 
evident that splitting the population before determination of the 
luminosity evolutions introduces a bias into the determinations, as 
discussed in ^4.21 

where j runs over all objects with a redshift lower than 
or equal to and m{j) is the number of objects with 
a redshift lower than the redshift of object j which are 
in object j's associated set. In this case, the associated 
set is again those objects with sufficient optical and ra- 
dio luminosity that they would be seen if they were at 
object j's redshift. The use of only the associated set for 
each object removes the biases introduced by the data 
truncation. Then the density evolution f(z) is 



fiz) = 



da{z) 
dz 



(13) 



(12) 



However, to determine the density evolution, the pre- 
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Fig. 7. — The cumulative density function a{z) vs. redshift for 
the RL (crosses) RQ (diamonds) and all (dashed line) quasars in 
the data set. The normalization of a{z) is arbitrary, and the RL 
data has been shifted vertically for clarity. A piecewise quadratic 
fit to a{z) is used to determine f(z) by equation 1131 



viously determined (in Q luminosity evolution must be 
taken out. Thus, the objects' optical and radio luminosi- 
ties, as well as the optical and radio luminosity limits 
for inclusion in the associated set for given redshifts, are 
scaled by taking out factors of ^opt(^) = (1 + z)^°p* and 
9rad{z) = (1 + ^j^''^^, with kopt and krad determined as 
above. 

Figures [7] and [8] show (7{z) and f{z) for the objects 
in the data set. We evaluate and display the density 
evolution separately for the RL and RQ objects and for 
the dataset as a whole to compare them. It is seen that 
the two groups, divided in this way, exhibit very similar 
density evolution. The number density of quasars seems 
to peak at between redshifts 1 and 1.5, a little earlier 
than general ly thought for the most luminous quasars 
(e.g. [Sh aver e t al. 1996), and earlier than that found in 
[Richards et al.» (^2QQ6i) , but similar to the peak f ound for 
less luminous q uasars by Hopkins ^et al.l (l2QQ7l ), and in 
agreement with Maloney & Petrosia^ (|1999 ^. 

Knowing both the luminosity evolutions gi{z)^ and the 
density evolution p(z), one can form the luminosity den- 
sity functions £i{z)^ which are the total rate of produc- 
tion of energy of quasars as a function of redshift. We 
show this for the radio luminosity density <£rad(^)- As 
evident the two populations of RL and RQ have very 
similarly shaped radio luminosity density functions (Fig- 
ure [9|) 

6. LOCAL LUMINOSITY FUNCTIONS 

6.1. General Considerations 

In a parallel procedure we can use the 'local' (redshift 
evolution taken out, or 'de-evolved') luminosity L[ distri- 
butions (and de-evolved luminosity thresholds) to deter- 
mine the 'local' LFs %jji{Li)^ where again the i represents 
the waveband, and the prime indicates that the luminos- 
ity evolution has been taken out. We first obtain the 
cumulative LF 
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Fig. 8. — The density evolution ^(z) vs. redshift for the for the 
RL (crosses) RQ (diamonds) and all (dashed line) quasars in the 
data set, shown with customary log scales. The normalization of 
fi^z) is arbitrary and the curves have been shifted vertically for 
clarity. 
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Fig. 9. — The radio luminosity density function i^rad(^) vs. red- 
shift for the for the RL (crosses) RQ (diamonds) and all (dashed 
line) quasars in the data set. The normalization of <£'rad(^) is arbi- 
trary, and the values have been shifted vertically for clarity. It is 
seen that the two populations have very similar luminosity density 
evolution with redshift. 



and, following iPetrosianI (|1992[ ). ^i{L'^ can be calculated 

by 



rik) 



) 



(15) 



where k runs over all objects with a luminosity greater 
than or equal to L^, and n(k) is the number of objects 
with a luminosity higher than the luminosity of object k 
which are in object A:'s associated set, determined in the 
same manner as above. The LF il^i{L[) is 



HL'i) 



dL' 



(16) 



In ^we have determined the luminosity evolution for 
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Fig. 10. — The cumulative local optical LF <l>opt(i^opt) 
quasars in the data set. A piecewise quadratic fit to ^L^^^^) is 
used to determine V^optC^opt) equation 1161 The normalization 
of <^»opt(^opt) ^^^^ arbitrary. 

two independent functions, the optical luminosity Lopt 
and the correlation reduced radio luminosity Lcrr- We 
can form the local optical ?/;opt(^opt) correlation re- 
duced radio V^crr(^crr) LFs straightforwardly, by taking 
the evolutions out. As before, the objects' luminosities, 
as well as the luminosity limits for inclusion in the as- 
sociated set for given redshifts, are scaled by taking out 
factors of ^crr(^) = (1 + ^)^'" and ^opt(^) = (1 + z^-^^ 
with kcrr and /Copt determined in ^ We use the notation 
L^L' = L/g{z). 

6.2. Local optical luminosity function 

Figures [TOl and [TTI show the local cumulative <l>opt(^opt) 
and differential tpopt{L^^ local optical LFs of the quasars 
in the I White et al.l (|2QQQf ) dataset, while figure [12] shows 
the local correlation reduced radio LF, ?/^crr(^crr)- 

The optical LF shows evidence of a break at 2 x 
10^^ erg see" ^JIz~^, winch was present already in data 
used in PetrosianI (|1973f ). Fitting a broken power law 
yields values — 2.0 i 0.2 and — 3.2 dz 0.2 below and above 
the break, respectively. If we allow for the possibility of 
additional uncertainty resulting from the consideration of 
possible radio incompleteness at faint fiuxes (see discus- 
sion in the range on the power law above the break 
increases to —2.8 ±0.4. As the optical LF has been stud- 
ied extensively in various AGN surveys, we can compare 
the slope of ?/^opt(^opt) obtained here to values reported 
in the literature. For example, 'Boyle et alj (|2000[ ). us- 
ing the 2dF optical data set (but with no radio overlap 
criteria) use a customary broken power law form for the 
LF, with values ranging from —1.39 to —3.95 for different 
realizations, showing reasonable agreement 

6.3. Local radio luminosity function 

With ?/^opt(^opt) '0crr(^crr)5 ^c cau determine the 
local radio LF '0rad(^rad) with equation [6l Figure [T3l 

It should be noted that they parameterize evolution differ- 
ently and work in absolute magnitudes rather than luminosities, 
however the slopes of their fits to the LF as they parameterize it 
are applicable, as can be seen in their section 3.2.2. 
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Fig. 11. — The local optical LF V^optC^opt) quasars in the 

data set. The normalization of '0opt(i^opt) here is arbitrary. 




Fig. 12. — The local correlation reduced radio LF '0crr(i^crr) 
the quasars in the data set. The normalization of ipcwiL^crr) here is 
arbitrary. For clarity, as plotted here, we have taken a numerical 
factor of (Lfid)°^ out of Lew (cf Equation O. 

shows the local radio LF '0rad(^rad) reconstructed in this 
way. It is seen that the local radio LF contains a possi- 
ble break around 10^^ ergsec"^ Hz~^, with a power law 
slope of —1.7 ±0.1 below the break and —2.4 ±0.1 above 
it. These ranges for the power law above the break are 
increased slightly to — 2.2 ± 0.3 if the effects of possible 
radio incompleteness are included, as in ^ The slope 
above the b r eak s een h ere is sim i lar to earlier results 
of Schmidt] (|1972r ) and iPetrosianI (|1973r ) which probed 
only those luminosities . A more complete comparison 
can be made with Mauch & Sadlerl (2007), who form ra- 
dio LFs of local sources in the Second Incremental Data 
Release of the 6 degree Field Galaxy Survey (6dFGS) ra- 
dio catalog. For the sources they identify as AGN, they 
find a break at 3.1 x 10^^ erg see" ^ Hz" ^, with slopes of 
-2.27±0.18 and -1.49±0.04 above and below the break 
(converting to luminosity units). 

7. DISTRIBUTION OF RADIO LOUDNESS RATIOS 

As stated in the introduction, naively one may expect 
that because the ratio R is independent of cosmological 
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Fig. 13. — The local radio LF V^radC^rad) quasars in the 

data set. The normalization of '0rad(^rad) ^^^^ arbitrary. 

model and nearly independent of redshift, the raw ob- 
served distribution would provide a good representation 
of the true distribution of this ratio. In Figure [14] we 
show this raw distribution by the triangles, arrived at by 
using the raw values of R from the data and forming a 
distribution in the manner of equations [15] and [16] with 
no data truncations. It appears that this naive approach 
shows a hint of possible bi-modality with \ogR ~ 1 as 
the dividing value 

As discussed in ^ we can reconstruct the local distri- 
bution of Gr'{R')^ as in Equation [8l which provides for 
a more proper accounting of the biases and truncations. 
The results of this calculation are also shown in Figure 
[TU The distribution calculated in this way clearly is dif- 
ferent than the raw distribution, and does not show any 
apparent bi-modality. There is still a possible feature 
in the same region {\ogR ~ 1) where the raw distribu- 
tion shows a dip. This feature is of marginal significance 
and results from the similarly shaped feature in ?/;crr(^crr) 
centered around L'^^^ = 10~^ergsec~^ Hz~^. Even if sig- 
nificant, this change in slope cannot be taken as evidence 
for two physically distinct populations, but could be a 
useful point to make an arbitrary division into RL and 
RQ objects. 

We also know that the best fit redshift evolution of the 
ratio, given equation[9l is gi^z) = (1 -\- z)"^-^. The change 
in the distributio n of with increasing redshift is also 
shown in Figure [T4F1 Another way to look at this is 
that we have found that the radio luminosity evolves at 
a different rate than the optical luminosity, with the con- 
sequence that their ratio is a function of redshift. The 
radio loudness of the population increases by a factor 
of 5 by redshift 1, and by a factor of 28 by redshift 
3. This is i n disagreement with the result presented by 
iJiang et "aTl ([2007) who show a decrease in fraction of 
RL sources with increasing redshift, which could be the 
case if the radio luminosity were to evolve more slowly 
than the optical luminosity. They however do not deter- 
mine individual evolutions or LFs. On the other hand. 

We note that in general apparent bi-modalities often do not 
stand up to rigorous statistical tests. 

Note that we have not included the density evolution which 
will shift the curves vertically but not change their shape. 
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Fig. 14. — The local distribution Gr{R) in the 5 GHz radio to 
2500 A optical luminosity ratio R, plotted as R x Gji(R), for the 
quasars in the data set. The stars are from Gji/{R') as determined 
by the method of Equation [S] taking account of the truncations 
and correlations in the luminosity evolutions, while the triangles 
result from forming a distribution with a naive use of the objects' 
raw ratio. The normalization is arbitrary and the curves have been 
shifted vertically for clarity. It is seen that the naive method gives 
a hint of a bi-modal distribution, while the proper method does 
not. Also shown is the proper radio loudness distribution Gji{R, z) 
at redshifts z=l (dashed line) and z=3 (dash-dot line), evolved 
according to the form of equation O 



iMiller et all dToOO") have noted that the fraction of RL 
quasars may increase with redshift, which they attribute 
to a difference in the evoluti o ns of the two populations 
(RL and RQ). Do noso et al. l (12009^ compute radio and 
optical LFs at differ ent re dshifts and reach the same con- 
clusion. Cirasuolo et alj (|2QQ6l) also find that the radio 
loud fraction may modestly incre ase at high redsh i ft. Al- 
though not directly comparable, iLaFranca et al.l ()2Q1Q| ) 
show a similar evolution for Rx, the ratio of radio to X- 
ray luminosity, as we show here for R. We note that our 
results favor one population, in the sense that the dis- 
tribution of G{R)^ recovered from considering the data 
truncations inherent in the survey and correlations be- 
tween the luminosities, is continuous. 

8. TESTS OF ASSUMPTIONS 

Power law parameterization: One may raise the con- 
cern that the simple power law parameterization used for 
the redshift luminosity evolutions (equation [5]) may not 
be the most ideal one. In particular, it may not accu- 
rately represent the evolutions at the highest redshifts 
considered here. To check this, we repeat the analysis 
with a different parameterization for the luminosity evo- 
lution which allows for a flattening at higher redshifts. 



hi 



(1 + ^) 



(17) 



where i again represents the optical or correlation re- 
duced radio luminosity. In this parameterization, the 
functional form for the radio luminosity evolution ^rad(^) 
and the evolution of the radio loudness parameter gR{z) 
are lengthier expressions involving /copt kcrr and a, given 
equations [3 and [9l 

This alternate parameterization for the evolutions does 
not appreciably effect the results. The best fit evolution 
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factors as a function of redshift under the alternate pa- 
rameterization differ very little from those in the simple 
parameterization to redshift 3.5 (the highest object in the 
sample). In a data set with higher redshift objects, the 
form of the parameterization will be more consequential. 

Luminosity dependent density evolution: Another con- 
cern may be that luminosity dependent density evolution 
(LDDE), which is not considered in the functional forms 
for the LF used here, may more accurately represent the 
evolution of the LF. As a check of this eflFect, we divide 
the data into high and low luminosity halves (cutting 
on optical luminosity), and check the similarity of the 
computed density evolutions for the two sets versus that 
computed assuming the absence of LDDE. Given that 
an artificial difference is already introduced in the two 
halves because there are a lack of low luminosity objects 
in the high redshift sample and a lack of high luminosity 
redshifts in the low redshift sample (see Figures [1] and 
[2]), we conclude that the density evolutions determined 
in this way are sufficiently similar to justify neglecting 
LDDE. 

Optical measurement errors: There is the possibility 
that errors in the optical magnitudes could lead to a bias 
which could affect the results. The bias introduced by 
these errors would be negligible if the LogN-LogS was ffat 
(i.e. N{> S) (X and dN/dS ex S'^). Since the num- 
ber density of sources increases with decreasing ffux, it is 
more likely that a source will be included than excluded. 
However, the magnitude of this effect will depend on the 
faint end source counts slope, and the shallower the slope 
the smaller the effect. The magnitude of this effect is 
proportional to [(1 -\-mbeiow) where S is the fractional 
error in flux and rribeiow is the faint end differential source 
counts slope (i.e. dN/dS cx S^^^^^"^). For the reported 
POSS-I magnitude errors of 0.2, for the faint end magni- 
tudes of ^17.8, 5 is less than 0.2. For miyeiow ^ —2 and 
^ ^ 0.2 the bias will be less than 4% averaged over the 
faintest fluxes, which would be manifest in raising the 
faint end cumulative source counts slope by an amount 
still considerably smaller than this , as it is fit over a 
larger range of fluxes. Previously, iCaditz fc PetrosianI 
(l993) investigated this effect for a flux limited dataset 
and showed that the difference was minor between using 
a Gaussian distribution of fluxes for each source instead 
of assuming a well deflned flux, for shallow faint end 
source counts such as here. 

Radio incompleteness: Lastly, the selection function for 
the FIRST objects in the White et al. sample used here 
might not be a sharp Heavyside function at a peak pixek 
flux density of 1 mJy, bu t rather smeared out. According 
to Figure 1 of iJiang eFal.l (^200 7). the selection function 
of FIRST for SDSS optically identifled quasars is such 
that at an integrated flux density of 1 mJy only about 
55% of sources are seen, and this number rises to 75% 
at 1.5 mJy and about 85% at 2 mJyEE This particular 
selection function would likely not be identical for the 
POSS-I optically identifled quasars of the sample we use 
here. Also, we have considered the sample to be limited 
by the peak pixel flux (i.e. surface brightness limited) 

The fuzzyness of the truncation boundary has a similar effect 
as the data measurement errors in the sense that it is unimportant 
for mijeiow = ~1 cind more important for larger deviations from 
this value. 



in the radio rather than being limited by the integrated 
flux, in accordance with the criteria set forward in White 
et al. So it is diflicult to directly compare the potential 
radio selection function here with the one in Jiang et al. 

The way to test the effects on our analysis is to repeat 
the analysis limiting the sample to a higher radio flux, 
where the sample would presumably be more complete, 
and determine the extent to which the calculated param- 
eters change in a systematic way. We have done so with 
lower radio flux limit of 2 mJy (486 objects), as opposed 
to the original 1 mJy (636 objects). The effect, propa- 
gated through the analysis, is primarily to extend the la 
uncertainties on /copt and ken in the direction of lower 
/Copt (1.25 on the low end) and higher /ccrr (3.75 on the 
high end), and to move the best bit values to 2 and 2.5 
respectively. There is no discernable effect on the value 
of the correlation parameter a. The modifled best flt val- 
ues and increased Icr uncertainty for A^opt 

and /Ccrr only 

slightly alters the la range and the best flt value of /crad, 
since the particular error ellipse shape means that lower 
/Copt values accompany higher /ccrr values. The main ef- 
fect on physical parameters then is to shift the best flt 
and low end la values for /copt downward. We also flnd 
that there is a negligible effect on the density evolutions, 
there is a small flattening effect on the high end opti- 
cal and radio LF slopes, which are reported in ^ and 
a negligible effect on the shape of the Gr{R) distribu- 
tion. To the extent that faint flux radio incompleteness 
is present in the sample considered here, it does not seem 
to have a large systematic effect on the determination of 
the parameters in this analysis. 

9. DISCUSSION 

We have used a general and robust method to de- 
termine the radio and optical luminosity evoluti ons si- 
multaneously for the quasars in the IWhite et aP ()2QQQl ) 
dataset, which combines 1.4 GHz radio and R-hand opti- 
cal data for 636 quasars ranging in redshifts from 0.02 to 
3.425 and over seven orders of magnitude in radio loud- 
ness. We flnd that the quasars exhibit more substantial 
radio evolution than optical evolution with redshift ( §4.21 
and Figure [6]). We also show that when divided into 
RL and RQ sets accordingly to the standard deflnition 
(divided by the value of the radio-loudness parameter 
R = 10), the two sub-populations exhibit similar den- 
sity evolution. The local optical and radio LFs that we 
obtain are consistent with previous determinations. 

Differences are noted with previous dete rminations of 
the ra dio luminosity evolution of quasars. IWillott et al.l 
()2001f ) also use a power law parameterization of the radio 
evolution with redshift for the radio bright sources they 
consider. Our result for the radio luminosity evolution, 
when evaluated for the data set as a whole is not con- 
sistent to within uncertainty w ith their results ( power 
laws ranging from 3.1 to 3.6). iStrazzullo erall (j2010f ) 
have recently obtained results with a radio survey to low 
(13.5 /iJy) flux limits, quoting k^ad = 2.7 ± 0.3 with the 
same parameterization used here for the sub-population 
that they identify as AGN. This is also not in agreement 
with the radio evolution determined in our analysis. In 
a future paper we intend to carry out the same analysis 
using the n i uch la rger SDSS sample of quasars used by 
'Jiang et al.J (j2007D . 

There has been much discussion as to whether RL and 



Quasar Luminosity Evolution 



11 



RQ quasars, defined solely by means of the radio loudness 
parameter as explained above, constitute a true contin- 
uum or two populations that can b e said t o be distinct 
in some way (e.g. [Ke ller man et alJ (| 19891) : llvezic et al.l 
(|2QQ2[ ): iCirasuolo et all (|2QQ3l) : llvezic e t al. (20041). Our 
analysis favors the former. First, we found that the di- 
vision of the quasar population into the two aforemen- 
tioned classes introduces strong biases into the simulta- 
neous determination of the radio and optical luminosity 
evolutions ( M.2\ see also Figure [6j). More importantly, 
as shown in ^ forming a distribution in the raw val- 
ues of the radio loudness parameter without taking 
into account the biases introduced by the truncations in 
the data and the correlated luminosity evolutions, results 
in a shape for the distribution which is very different 
from the true distribution; e.g. it shows a possible dip in 
Gr{R) dit R ^ 10, while the true distribution is rather 
smooth and shows at most a modest feature. Even if 
the feature is real, it does not suggest any bi-modality, 
but rather a continuous range of physical properties in a 
single population. We also find that at higher redshifts 
and optical luminosities the radio loudness become more 
pronoun ced. This is opposite of the trends presented by 
Jiang et a l. (2007), but in agreement with others such as 
Donoso et al.i (;2009i ). 

Accessing the unbiased distribution of the radio loud- 
ness parameter for quasar sources is crucial not only for 
understanding the cosmological evolution of this class of 
active galaxies, but also for understanding jet launching 
processes in the vicinities of supermassive black holes. In 
this context, we note that the observed optical fluxes of 
quasars are dominated by the emission from an accretion 
disk accreting at relatively high rates, around 1% — 100% 
of Eddington, and therefore radiating with 10% effi- 
ciency. Hence, the optical luminosity is a very good mea- 
sure of the total accretion power in quasar sources. On 
the other hand, the observed radio fluxes of the discussed 
class of objects are expected to originate in the outflow- 
ing magnetized plasma. In particular, the radio emis- 
sion of quasars is produced predominantly via the syn- 
chrotron emission of relativistic well-collimated jets (in 
the case of very radio loud sources), or via the cyclotron 
and/or free- free emission of at most mildly-relativistic 
disk winds (in the case of very radio quiet nuclei). In 
both cases, the observed radio luminosities should be 
considered as proxies for the kinetic luminosities of the 
outflowing matter. Therefore, the radio loudness R char- 
acterizes the efficiency of the production of jets/outflows 
for a given accretion power. 

The lack of any clear bi-modality centered around 
R=10 in the distribution of the radio loudness parameter 
for quasars, as advocated here, implies then that there is 
no critical change in the parameters of the central engine 
between the RL quasars (those producing extremely pow- 
erful relativistic jets), and the RQ ones (those produc- 
ing only mildly-relativistic and uncollimated disk winds). 
This is a crucial piece of evidence for understanding still 
debated mechanisms for jet launching in AGNs. Note, 
for example, that our flnding is hardly consistent with 
the idea that RL quasars possess counter-rotating (with 
respect to the black hole spin) accretion disks, as opposed 
to RQ quasars with co-rotating disks only. Instead, the 
nuclei of jetted and non-jetted quasar sources seem to 
be intrinsically very similar, differing only smoothly and 



continuously in some particular respects. But it has to be 
emphasized that here we do not discuss the whole pop- 
ulation of AGNs (including, e.g., Seyfert galaxies) but 
only strictly the quasar populati on (see the related recent 
discussion in Sikora et al. 2007", and references therein). 
We note that our analysis does not address the question 
of whether there may be a large population of quasars 
with values of R beyond the range of the present sample, 
in particular we cannot rule out a signiflcant population 
with very low values of R {logR < 0.01). 

Another (possibly related) result we flnd is that the 
radio loudness may increase with increasing optical and 
radio luminosities, as the best flt value for the correlation 
parameter a is 1.3, although a strictly linear correlation 
is not ruled out with much signiflcance, and furthermore, 
one may dispute that the correlation is inherent in the 
population (see the discussion in Appendix B). However, 
if this super-linear correlation were indeed intrinsically 
the case, it would imply the existence of some connec- 
tion between the efficiency of formation of relativistic 
jets and accretion power, which may in turn depend on 
the combination of the evolving accretion rate and black 
hole spin (see in this context Tchekhovskoy et al. 201Q|), 
given the observed increase in radio loudness with red- 
shift. Note, however, that even though we have used 
simple one-parameter functions to describe the emerg- 
ing correlations, it is possible that some of them are 
more complex. For example, the correlation index a be- 
tween the radio and optical luminosities may be close to 
unity only for low luminosity objects, but much larger 
than that for more luminous (and therefore radio loud) 
quasars. More data and further analysis is needed to ad- 
dress this and similar issues, which may provide further 
constraints on theoretical models. 

Another application of the presented analysis is re- 
lated to the understanding of the origin of the cosmic 
background r adiati on in the radio frequency regime. In 
iSingal et al.l (|20Toh we estimated the fractional contri- 
bution of quasars to the Cosmic Ra dio Backgroun d, as- 
suming the level reported by Fixse n et al.l (l2010f ). In 
general the flux of the objects fltting the deflnition of 
radio loud is well characterized by current interferomet- 
ric radio surveys so that their contribution to the total 
radio background intensity can be estimated to be 15% 
to 25% of the observed value. In the earlier work we 
also estimated the total contribution to the background 
of the radio quiet objects, and found it to be between 
1% and 2% for favored models of quasar luminosity evo- 
lution. This estimate was based on integrating values 
of the quasar bolometric LF, as reported in the litera- 
ture, over redshift, applying a mapping between optical 
and radio luminosity, and assuming that the optical and 
radio luminosities had identical redshift evolutions. We 
also noted there that the contribution we estimated was 
dependent on the later assumption and would be revised 
in the case of differing radio and optical luminosity evo- 
lutions. As we see here that the population of quasars 
has greater radio evolution relative to optical, the con- 
tribution of RQ quasars to the radio background will be 
somewhat larger than the value reported in our previous 
work. We will present a quantitative determination in a 
forthcoming paper. 
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APPENDIX 
APPENDIX A 

Distribution of Radio Loudness R : As described in the introduction one of the aims of this paper is deter- 
mination of the distribution of the radio to optical luminosity ratio R = i^rad/^opt from an observed sample of radio 
and optical fluxes and redshifts. This requires a proper accounting for correlations and evolutions of optical and 
radio luminosities and the observational selection effects. Here we describe the how these factors affect the observed 
distribution of R. 

The true or intrinsic distribution of R values is related to the radio and optical LFs ^(Lopt, ^rad, ^) as 

Gt{R,z) = j "^{Lopt.R Lopt,z)LoptdLopt = J ^ i^rad, ^rad (1) 

(compare to the separable, local form in Equation [8]) 

The observed distribution on the other hand is different because the observational selection effects truncate the 
data. For example for a sample with weh defined flux limits /rad > fm,rad, /opt > fm,opt, and i^obs = fm,rad/ fm,opt 
the observed distributions is 

POO 

Gohs{R, Z)= j ^(i^opt, R ^opt, Z) Lopt dLo^t for R > Robs (2) 

GobsiR, Z)= I * f Lrad, z\ 



I ^^,^rad,^ ) ^rad 

.,rad(2) 
pOO 

= / ^(^opt, R ^opt, z) Lopt dLopt for R < i?obs, (3) 

where 

R = (i^opt/^rad)^, ^min,opt(^) = d\{z) K o^t f m ,opt and Lmin,rad(^) = ^7rdl{z)Kradfm,rad' (4) 

and Ki is the K-correction factor for waveband i. All of these are obtained from the observed distributions of the 
fluxes and redshifts. If we approximate the observed distribution of fluxes and redshifts by a continuous function 

^obs(/opt,/rad,^) then 

-^oDt -^^rad \ f ^ \ 1 



*(Lopt, L,ad, z) - nobs ( ^^al{z)Ko,, ' 47r4(^)ifrad ' V Wdl ) Ko,,K,,,V'{z) 



SO that 



GT{R,z) = {KoptKrs.dV')~^ I nobs(/opt,^/opt,^)/opt4fopt (6) 

JO 

POO 

Gohs{R,z) = {KoptKrs.dV')~^ I nobs(/opt,^/opt,^)/opt4fopt, (7) 

where V = dV{z)/dz, and /^m = fm,opt for R > Robs and /^m = fm,opt{Rohs/R) for R < Robs- Note that R depends 
on redshift to the extent that the optical and radio K-corrections are different. However, since the radio and optical 
spectra can be well approximated by power laws with the almost same spectral index the two K-corrections are almost 
equal. In what follows we ignore this small difference and set R = R. The sources with R > Robs can be called 
optically limited because their optical flux is the main determining factor for their inclusion in the sample. Similarly 
sources with R < i^obs can be called radio limited. 

Clearly the intrinsic and observed distributions of R are different. In reality the situation is slightly more complicated 
because the lower limit of the integration does not extend to zero. The samples of available quasars are truncated 
also by minimum luminosities, say Lm,opt and Lm,rad which introduces a second critical value for i?, namely i^int = 
Lm,rad/Lm,opt' The abovc equations are valid for redsifts z > 2;niin,opt or 2;min,rad defined as 

^min,i)-^i(^min,i)/m,i- (8) 

For z < Zmin (defined as the lower of 2;min,opt and ^min,rad) there is no truncation due to fiux limits and Gobs = Gt- 
It is convenient to define 

POO 

-^opt? ^) -^opt dL^pt (9) 

J i 

SO that the true distribution can be written as 

GT{R,z) = ^R,z;Lm,opt) for R>Rint (10) 
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Gt{R,z) = '^>(^R,z;^^^!^^^^ for R < R^^f (11) 

The observed distribution depends on the relative values of i^int and i^obs- 

For a sample where Robs > ^int, which means that 2;min,rad > ^min,opt, such a sample may be classified as mainly 
optically selected. (For purely optically selected sample i^obs = 0-) In this case observed distribution will be the same 
as the true distribution (Gobs = Gt) for all redshifts z < ^min,rad- Obviously, the reverse is true in the opposite case 
with radio exchanged for optical. At higher redsifts 

GobsjR^z) ^{R,z;Lmin,opt{z)) 

7^ ~ — ^TTT^ ? ^ — ^ > ^int U^j 

/TS / D ^. -Robs ^min,opt(^) ^ 

u } ^^^^^^^ ^^^1 



Gt (h ( T) 'V' -Rint Lm,opt 

R 



A simple Example: Let us assume that the radio and optical luminosities are uncorrelated and do not evolve so 
that we have ^(Lopt, i^rad, ^) = '0opt(^opt) '^rad(^rad) p{z)^ whcrc p{z) describes the density evolution. Furthermore if 
we assume simple power law LFs '0opt(^opt) = ^opt^op^°'^* '0rad(^rad) = ^rad^^^d'^'''^, it is easy to show that 

GT{R,z)(x{R/Rint)^-''^'''' for R>Rint and Gt{R, z) oc {R/ RintT°'''~^ for R < Rint. (14) 
Similarly, it is easy to show that for redshifts z > Zmin 

/ T (^\\^ f 1 for > Rohs > Rint, 

GoM z)=Gt[ ) X \ {Robs/R)^ for Robs <R< RinU 

^ ^^'^^^ ^ [{Rohs/Rint)^ for Rol,s < Rint < R 



where /3 Ofopt + Q^rad - 2. 

The observed and true distributions have different shapes in the range between the intrinsically and observationally 
limiting values of R. The shapes become identical when these two values are equal. Larger differences will be the case 
if the radio and optical luminosities are correlated non-linearly and undergo different kinds of luminosity evolution. 
As a result the fraction of RL or RQ sources (arbitrarily chosen at some value or R) will vary with redshift and/or 
luminosities. For example if the radio and optical luminosities were correlated linearly making the the radio loudness 
independent of both luminosities then we will be dealing with a separable LF ?/^(i?, Lopt, ^) = G {R)iIj opt {L opt ^ z). In 
such a case, because of fiux limits at any redshift the observed range of extends to i^min(^opt7 z) which decreases with 
increasing Lopt but increases with increasing z. This will cause the fraction of RL sources to decrease with luminosity 
but increase with redshift. All such trends can be determined by proper accounting of the correlations and evolution 
as described in this paper. 

APPENDIX B 

Lum inosity correlations and flux limits: As pointed out by the referee, it has been suggested by IAntonuc"cj| 
()2Qlll) that the observed correlation betwee n the luminosities may b e induced by observ ational sel ection effe cts. Other 
works that have examined the issue include [Khembavi et al.l (1 19861 ), iFeigelson fc BergI ([T983) and iChanan I ([1983). It 
is not clear whether the correlation seen in this work between Lrad and Lopt is inherent in the quasar population, 
or is introduced by the selection effects of the surveys. In some sense this question does not matter for the analysis 
presented here, because the rotation to Lcrr is a technique required to achieve independent variables {Lopt and Lcrr) 
in the context of the data present so that we can recover the inherent redshift evolutions. This is independent of 
the underlying luminosity-luminosity correlation of the population. However, it is of general interest whether the 
correlation seen between radio and optical luminosity is inherent. 

In order to begin an investigation of this, we simulate via Monte Carlo (MC) techniques two cases of quasar 
populations, distributing the objects in radio luminosity, optical luminosity, and redshift, and applying the fiux limits 
of the White et al. dataset to achieve for both MC populations 'observed' datasets of a similar size to the real White 
et al. data used in this analysis. For the MC populations we use a redshift distribution with f(z) ex to redshift 1.2 
and then constant density with redshift for 1.2 < 2; < 3.2. In the first, uncorrelated, case, we draw the optical and 
radio luminosities of the population from separate distributions with power law slopes in luminosity of —2. In the 
second, correlated, case, we draw the optical luminosities of the population from a distribution with a power law slope 
in luminosity of —2, then the radio luminosities of the population are assigned according to Lrad ^ ^opt^'^ and then 
randomized about that value by a factor of 10^-^^^ where A is a normally distributed random number with mean of 
zero and standard deviation of one. For this simple analysis we have not included the effects of luminosity evolution, 
because it would require an orders of magnitude larger simulation. 
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Fig. 15. — Raw observed cumulative distribution in redshift of the White et al. dataset (small plusses), the 'observed' dataset from 
the uncorrelated MC population (diamonds), and the 'observed' dataset from the correlated MC population (stars). The normalization is 
arbitrary. The redshift distribution of the observed correlated MC population much more closely resembles that of the real observed data. 

Our results show that that indeed the radio and opti cal lu minosities of the 'observed' datasets in both cases are 
correlated. An analysis identical to that presented in ^4.11 reveals that the observed data for the correlated MC 
population has a value for the correlation index a of 1.3 ± 0.3 while that of the uncorrelated MC population has a 
value of 0.85 ± 0.15. However, these 'observed' sets can be compared to the real White et. al dataset, and it is seen 
that the observed dataset for the correlated MC population more closely resembles the real dataset than that of the 
uncorrelated MC population does, for instance comparing the redshift distributions of the observed objects, shown in 
Figure [m From this analysis we conclude that the correlation seen in the White et al. dataset between Lrad and Lopt 
could be inherent in the population. Exact determinations of what fraction of the observed correlations is inherent in 
the population and what fraction is due to selection effects depends on the values of the many parameters describing 
the LFs and evolutions. This is an important issue in many areas of astrophysics and requires considerable work which 
is beyond the scope of this paper. In a forthcoming work we will address this general question. 



